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Abstract 

Discussed in the paper is a mixed scenario of the charged liquid surface re- 
construction in the situation where the 2D surface charge density is close to its 
saturation value. 

One of the instabilities studied in classical hydrodynamics is the Frenkel-Tonks (FT) 
instability [1,2] arising in a threshold manner on the charged liquid surface and resulting in 
its deformation. In contrast to other known instabilities, such as the Rayleigh instability 
of a cylindrical jet [3], von Karman trace generated by moving cylinder (sphere) [4], Taylor 
vortex instability of a viscous layer between two co-axial rotating cylinders [5], etc., the 
FT decomposition process can be halted a new metastable state with finite corrugation 
depth can be formed (charged liquid surface reconstruction). The surface reconstruction 
was studied by different authors at different times [6-19]. However, by the present time a 
consistent comprehensive picture of the instability development is only available for small 
values of the parameter v <C 1 

u = n s /nT ax , nr x = V^P9, (1) 
(p and a are the liquid density and surface tension, g is the gravity acceleration, k = 



pg/a is the capillary length), determining the liquid boundary population by localized 
charges. Explicitly manifested here is the dual nature of the FT instability development: 
the spectral scenario where the frequency imaginary part changes its sign as a function of 
external parameters or, alternatively, or the surface corrugation arises due to fluctuations. 
Under these conditions it is possible to draw the phase diagram branches in the (z/, EJ) 
plane operating with the concepts used in the physics of the 1-st order phase transitions: 
spinodal, binodal, etc. [20]. Here E- is the electric field above the charged liquid surface. 
At the v — 1 the electric field is E_ = 0. 

For finite values v < 1 available in the theory are only isolated facts which do not form 
any coherent picture. For example, it is easy to extend the the spinodal in the (v,EJ)- 
plane into the v < 1 domain. Here the alternative part of the problem, the nucleation- 
controlled "decomposition" of a flat charged surface (the binodal) is left in the uncertain 
state (there are no typical scale for the charge involved into formation of a dimple). On 
the other hand, the experiment unequivocally indicates a periodic reconstruction [18,19]. 
This difficulty is lifted if one considers the FT instability temporal evolution. Authors 
of Refs. [21,22] noted that the spinodal scenario of the surface decomposition with a 
fixed total charge leads to the appearance of a periodic array of neutral spots on charged 
surface. This behavior obviously means the end of spino- and the beginning of binodal 
evolution with well defined initial conditions: the charge per arising dimple is known 
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and all the dimples are identical. We refer to the scenario beginning with the "spino- 
" regime and continuing with the "bino-" evolution of the charged surface deformation 
as the mixed one. This scenario reproduces in a qualitatively correct way basic known 
properties of the FT instability in the v < 1 domain: periodicity of the reconstruction, 
loss of equipotentiality along the corrugated liquid surface, and jump-like behavior of the 
corrugation amplitude when the stationary reconstructed liquid boundary is formed. 

One should emphasize that the mixed decomposition scenario of two-component sys- 
tems is a rule rather than exception, the point that has not yet been fully taken advantage 
of. Consider, for example, the phase diagram in the vicinity of its extremum for a two- 
component solution [23]. Standard system transition from a "stable" point a to unstable 
point b by a jump of temperature initiates the onset of the spinodal decomposition in the 
solution along the c — b — c\ isotherm. In this process the system undergoing decompo- 
sition in the "directions" from b to c, C\ approaches the spinodal-binodal boundary from 
the inside of the spinodal zone rather than from the uniform state which takes place in 
the standard transitions of the d — e type. Starting from the c* — c* points the evolution 
becomes binodal one characterized by the spatially modulated stratified density arising 
at the end of the spinodal decomposition phase (just like in the above reconstruction 
problem) rather than by a random nucleation-controlled density. 

Turning back to the charged liquid surface, let us formulate a number of arising prob- 
lems. First of all, it is highly desirable to make sure that the stationary "dimpled" state 
in energetically more favorable than the unstable flat state (experimentally, this is quite 
obvious since no deformation would otherwise be observed), one should also estimate the 
charge localization degree in the sense of the inequality 

kR < 1, (2) 

where R is the single charged spot radius. 

1. In the standard flat capacitor geometry with the distance h between the metal 
electrodes and the liquid film of thickness d covered by a charged 2D system characterized 
by the number density n s the problem of finding the surface properties is assumed to be 
"bulk-like" if h > d » a, where a = k" 1 is the so-called capillary length. In addition, 
we are only interested in the saturation case where the electric field E_ above the liquid 
is zero while the density n s has its maximum possible value, i.e. v — 1. Under these 
conditions the electric field E™ ax in the liquid film also reached its maximum value equal 
to 

E\ ax = 4iren™ ax (3) 

where v is taken from (1). 

It is convenient to compare the system energy in different states (flat W and corrugated 
W) by first considering the collective effects and then calculating the contribution of a 
single Wigner-Seitz cell. There are two collective effects. The first one is the gain in 
Coulomb energy Aw c due to charge ordering. This contribution plays the major role in 
Coulomb crystallization on flat substrates and is estimated to be [24] 

Aw c ~ -l.lQ 2 /a (4) 
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where Q is the single dimple charge and a is the lattice spacing. 

The second effect is specific to to the charged liquid surface. A continuous electron disk 
of radius pressed by the electric field E™ x to the liquid surface squeezes some liquid into 
the neutral neighborhood of its perimeter. This effect raises the net electron system energy 
by AW^ (see Appendix). The same disk driven to the lattice state displaces a different 
amount of liquid (actually, it will be shown below that in that case no liquid is displaced 
it all). Therefore, the surplus energy AW% favors the charged surface reconstruction. 

In terms of Aw c and AW^ the general requirement 

W > W 



is written as 

w c + AW% > w c + nwi, w c — w c = nAw c , AW^ — nAw^. (5) 

Here w c is the total Coulomb energy of a capacitor with a flat charged liquid layer, w c is the 
same energy in the presence of a cluster lattice, wi is the total energy of a single multiply 
charged dimple, n > is the total number of such dimples arising in the reconstruction 
process. When writing Eq. (5) we have assumed that in the mixed reconstruction the 
corrugation possesses the following property: each cell of the arising structure contains a 
single multiply charged dimple with a charged core of radius R such that a/R 1 and 
the array of dimples has a spatial period comparable to the capillary length. This means 
that the inequality (5) we are interested in is satisfied if 

wi < Aw c + Avj£. (6) 

To estimate wi one could employ the available results [22] on the properties of indi- 
vidual dimples. In that case the energy wi becomes zero at 

E crit = l.UQE™ x 

which is qualitatively quite acceptable. However, the product kR in that case yields 
kR ~ 1.5 > 1, making the single dimple approximation completely meaningless in the 
calculations of w\ for mixed reconstruction. Here the problem should be reformulated 
from the very beginning taking into account the finite size of the cell containing a single 
charged spot. 

The appropriate approach is given by the Wigner-Seitz (WS) approximation with the 
external cell radius of the order of a and the charged spot radius (actually, there exist two 
alternative distributions differing by whether charged or empty spots occupy the central 
part of the WS cell; qualitative estimates given in the Appendix favor the charged central 
spots). In that case the following set of equations is to be solved: 

^ x 9 f n(r')r'dr' 

eE ± Z(r) + 27re 2 / ; = A e (7) 
J \r — r I 

- K 2 i - P el {r)/a = \, P el {r) = eE ± n(r) (8) 



3 



d£/dr\ r ^ a = 0, 2tt rdr£(r) = (9) 
Jo 

2tt J n(r)rdr = N, N = ira 2 n s (10) 



Here £(r) is the self-consistent helium surface deformation, k = ypg/a is the liquid 
capillary length, P e i{r) = eE±n(r) is the electron pressure on the helium surface, E± is 
the electric field pressing charges towards the liquid surface (in our case, according to Eq. 
(3), E± = E + = 4iren™ ax ), A e is the Lagrange multiplier accounting for the fixed total 
charge (10). Physically, this multiplier is equivalent to wi (6). A^ is another Lagrange 
multiplier ensuring the total liquid volume conservation (9). 

Assuming the electron spot radius R to be sufficiently small, (kR < 1), it is convenient 
to represent the surface deformation £(r) in the central part of the cell as a series 

£(r)^(0) + ^"(0)r- 2 + ... (11) 

In that case, by analogy with the contact Hertz problem [25] 

E ^ = l Ne fw^ (13) 

X e -eE ±m = - 4 N e >l w - j7 , (14) 

In addition, in the central part of the cell the Eq. (8) becomes 

f"(0) - /t 2 £(0) - A f = 3QE ± /(4irR 2 a), (15) 

Relations (13) and (15) contain there dimple characteristics: £"(0), R, £(0). To close 
this set of equations, one should add a definition of £(0), i.e. the solution of (8) which can 
be obtained by expanding £(r) into a series in appropriate Bessel functions. By writing 

f(r) =£ + f(r), £ = -A € , (16) 

and choosing the functions 

dJ {fi n r/a) 
J {n n rj a) , — 1 r=a = 

as a basis set so that to automatically satisfy the first boundary condition in (9) one can 
obtain to solution to Eq. (8) as a series 

oo 

e(r) = e + E& J o(Mnr/a), (17) 

T( = l 
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~ Af(/i n ,R/a) , . sinx — xcosx A 2>QE± 

^ n = ~ 1 .9 tTT. V' J\ x ) 



The quantity £ here is obtained from the second boundary condition (9) and proves to be 
zero. 

Analysis of the equation set (13), (15), (17) allows one to suggest (see arguments cited 
in the Appendix) that at the initial stage of the dimple reconstruction development the 
quantity £ (0) does not play any significant role in these equations. Hence, the relationships 
between £"(0) and R immediately yield a closed equation set from which they can be easily 
found: 

^(0) = f^jp R = * 2 */Ei, (18) 

Here the quantity £(0) is obtained from Eq. (17) with R calculated from Eq. (18). 

The final stage in the estimation of the mixed scenario feasibility is the estimation of 
parameters which are most critical for its realization. By setting in the general formulas 
(7)-(18) E ± = E + = 4Tien™ ax , one has 



0.4 < 1, scmv ' ~ 0.3 < 1, \T lt = w? nt < A<*' + A< (19) 
a a 

Here n™ ax is taken from (1) and w\ from (6). 

Estimates (19) reveal that the mixed scenario of the charged liquid helium surface 
reconstruction seems to be rather acceptable in the vicinity of the filling factors (1) close 
to unity. In the present paper, it is impossible to obtain more definite conclusions since the 
employed approximation (B-3) cannot claim any quantitatively correct statements. We 
do not see any other approaches (different from the mixed scenario) which could proved at 
least qualitative explanation of the observed periodic reconstruction of the liquid helium 
surface in the range of v < 1. 

This work was partly supported by the RFBR. 
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Appendix 

1. Let R* be the charged area radius and L > the radius of entire liquid surface 
between vertical walls. The equilibrium of the cell as a whole is defined by the following 
equations: 

pgZo + Pei = pgti, (Pi) 
R% + (L 2 -Rl)ti 1 = 0, L > R*. (P2) 

Here £ is the liquid surface deformation in the charged area, a, £i is the the liquid surface 
deformation outside the electron disk, and P e \ is the effective pressure in the charged area. 
Finding from (P2) the quantity £i 

6 = -^Rl/{L 2 - Rl) 
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and substituting it into (PI) one has 



Z = -P a /g m , 9* = g(i + T ^ ir2 ) (P3) 

The additional energy AW^ due to deformation £ has the scale 

AW 6 ~ nRlP^o (PA) 

Combining all the necessary definitions one obtains from (P4) 

Anna 2 R 2 
AW^nA W( , A w , = 7 ^- Iy , m ,L) = (1 + w ^). (PA) 

This expression is used for estimates in the main text of the paper. 

2. Within the alternative including either charged or empty spots, it is convenient to 
start from qualitative estimates. Suppose we consider a cylindrical Wigner-Seitz model 
with the spot external radius of the order of a, charged spot radius R < a, and charge- 
free ring with area 7r(a 2 — R 2 ). If the last formula is interpreted as an integral relation, 
one can estimate within the cylindrical model the probability of scenario described in 
Ref. [21] where the initial stage of the reconstruction is assumed to be represented by a 
periodic array of empty (charge-free) spots whose area is different from zero to the extent 
of positive supercriticality SE = (E_ — El), 

5E=(E„-E C L). (P5) 

However, this is most probably not so because of the following reasons. 

Within the WS cell both the local equilibrium conditions and liquid volume conserva- 
tion should be satisfied, just as in (PI). Then, by similar arguments one has 

io = -Pjg\ g* = g(l + ^-=) (P6) 

cr — R z 

According to (P6), development of a stationary corrugation under the conditions n(a 2 — 
R 2 ) — > (empty spots with low area) is unlikely since in that domain the effective value 
of g* — > oo. In that case the appearance of corrugation, i.e. growth of £ (P6) is hindered. 
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